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Abstract 

A deformation of the classical trigonometric BC,i Sutherland system is derived via 
Hamiltonian reduction of the Heisenberg double of SU(2n). We apply a natural Poisson- 
Lie analogue of the Kazhdan-Kostant-Sternberg type reduction of the free particle on 
SU(2re) that leads to the BC„ Sutherland system. We prove that this yields a Liouville 
integrable Hamiltonian system and construct a globally valid model of the smooth re¬ 
duced phase space wherein the commuting flows are complete. We point out that the 
reduced system, which contains 3 independent coupling constants besides the deforma¬ 
tion parameter, can be recovered (at least on a dense submanifold) as a singular limit of 
the standard 5-coupling deformation due to van Diejen. Our findings complement and 
further develop those obtained recently by Marshall on the hyperbolic case by reduction 
of the Heisenberg double of SU(n,n). 
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1 Introduction 

Models amenable to exact treatment provide key paradigms for our understanding of natural 
phenomena and form a fertile field of research crossing the border of physics and mathematics. 
The study of integrable Hamiltonian systems is a very active subheld with particularly strong 
ties to group theory and symplectic geometry. For reviews, see e.g. [nii22i[3niEiis]- One of the 
time-honoured approaches to such systems consists in viewing them as ‘shadows’ of natural 
free systems enjoying high symmetries. This is alternatively known as the projection method 
or as Hamiltonian reduction 121ES]. The list of the free ‘master systems’ is monotonically 
expanding in time. To name a few, it includes free particles on Lie groups together with their 
Poisson-Lie symmetric deformations and quasi-Hamiltonian analogues. For example, it was 
shown in the pioneering paper d that the integrable many-body system of Sutherland [31] , 
which describes particles on the circle interacting via a pair potential given by the inverse 
square of the chord-distance, is a reduction of the free particle on the unitary group U(n). 
Various deformations of the Sutherland system due to Ruijsenaars and Schneider [H [29| 
were derived mm from Poisson-Lie symmetric free motion on U(n), whose phase space 
is the Heisenberg double [33] of the Poisson-Lie group U(n), and from the internally fused 
quasi-Hamiltonian double [2] of U(?7,), which arose from Chern-Simons held theory. 

The projection method was enriched by an interesting recent contribution of Marshall 
poi. who obtained an integrable Ruijsenaars-Schneider (RS) type system by reducing the 
Heisenberg double of SU(n, u), which directly motivated our present worl0. Here, we shall 
deal with a reduction of the Heisenberg double of SU(2?7,) and derive a Liouville integrable 
Hamiltonian system related to Marshall’s one in a way similar to the connection between the 
original trigonometric Sutherland system and its hyperbolic variant. Although this is essen¬ 
tially analytic continuation, it should be noted that the resulting systems are qualitatively 
different in their dynamical characteristics and global features. In addition, what we hope 
makes our work worthwhile is that our treatment is different from the one in [23] in several 
respects and we go considerably further regarding the global characterization of the reduced 
phase space and the completeness of the relevant Hamiltonian flows. 

The main Hamiltonian of the system that we obtain can be displayed as follows 


H{p,q-,x,u,v) = -^ 


e-2- + e2- " 


g-2« ^ ^2v 



k=l 


Here u, v and x are real coupling parameters that will be assumed to satisfy 

u < V, V ^ —u and x ^ 0. 

The components of q parametrize the torus T„ by and p belongs to the domain 
Cx := {p e M” \ 0 > Pi, Pk - Pk+i > \x\/2 {k = 1,... ,n - 1)}. 

The dynamics is then dehned via the symplectic form 


( 1 . 2 ) 


(1.3) 


n 



^The relation is ‘symmetric’ as the problem studied by Marshall was originally suggested by one of us. 
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It will be shown that this system results by restricting a reduced free system on a dense open 
submanifold of the pertinent reduced phase space. The Hamiltonian flow is complete on the 
full reduced phase space, but it can leave the submanifold parametrized by Cx x T„. By 
glancing at the form of the Hamiltonian, one may say that it represents an RS type system 
coupled to external helds. Since differences of the ‘position variables’ pk appear, one feels 
that this Hamiltonian somehow corresponds to an A-type root system. 

To better understand the nature of this model, let us now introduce new Darboux variables 
Qk, Pk following essentially pni as 

exp(pfc) = sin(gfc) and = Pfc tan(gfc). (1.5) 

In terms of these variables H{p, q; x, u, v) = x, u, v) with the ‘new Hamiltonian’ 


ni{q,p;x,u,v) = 


e-2“ + e2" - 


V- 

Cl 


^ sin (qj) 


5^cos(pjtan(gj)) 

i=i 


1 + 

1-• 2( \ -^ 

sm [qj) 


^^2{v-u) 


X 


n 


k=l 

(Mi) 


4sin^(gj) — sin^(2gj) 

2sinh^ (l) sin^(gj) sin^(gfc) 
sm^((jj - + qt) 


( 1 . 6 ) 


Remarkably, only such combinations of the new ‘position variables’ qk appear that are natu¬ 
rally associated with the BC„ root system and the Hamiltonian "Hi enjoys symmetry under 
the corresponding Weyl group. Thus now one may wish to attach the Hamiltonian "Hi to the 
BC„ root system. Indeed, this interpretation is preferable for the following reason. Introduce 
the scale parameter (corresponding to the inverse of the velocity of light in the original RS 
system) /? > 0 and make the substitutions 


u —>• f3u, V —>• /3v, X —)■ f3x, p —)■ f3p, u —)■ pCj. 

Then consider the deformed Hamiltonian 

ni 3 {q,p;x,u,v) := 'Hi{q,/3p]/3x,/3u,/3v). 

The point is that one can then verify the following relation: 

Hp{q,p-,x,u,v) -n snthf ^ 

^ P; 7,7i, 72 ), 


where 




j = l 


7 


7 


sm\qj-qk) siuRg^ + g^) 


(1.7) 


( 1 . 8 ) 


(1.9) 


^sin"(2gj) 


72 


( 1 . 10 ) 


is the standard trigonometric BCn Sutherland Hamiltonian with coupling constants 


X 


7 = ^, 7i = 2nn, 72 = 2 {v - uf. 


( 1 . 11 ) 
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Note that the domain of the variables g,p, and correspondingly that of q, p also depends on 
/3, and in the /3 —?• 0 limit it is easily seen that we recover the usual BC„ domain 

qi> q 2 > ■■■> qn>0, p G M”. (1.12) 

In conclusion, we see that H in its equivalent form l-ip is a 1-parameter deformation of the 
trigonometric BC„ Sutherland Hamiltonian. We remark in passing that the conditions fll. 2 p 
imply that 72 > 0 and dyi -|- 72 > 0 , which guarantee that the flows of complete 

on the domain fll.l 2 p . 

Marshall ra obtained similar results for an analogous deformation of the hyperbolic BC„ 
Sutherland Hamiltonian. His deformed Hamiltonian differs from fll.ll) above in some impor¬ 
tant signs and in the relevant domain of the ‘position variables’ p. Although in our impression 
the completeness of the reduced Hamiltonian flows was not treated in a satisfactory way in 
[2U] . the completeness proof that we shall present can be adapted to Marshall’s case as well. 

It is natural to ask how the system studied in the present paper (and its cousin in [20]) is 
related to van Diejen’s [35] 5-coupling trigonometric BC„ system? It was shown already in 
[35] that the 5-coupling trigonometric system is a deformation of the BC„ Sutherland system, 
and later [SB] several other integrable systems were also derived as its (‘Inozemtsev type ’ [IB]) 
limits. Motivated by this, we can show that the Hamiltonian fll.ip is a singulai[§ limit of van 
Diejen’s general Hamiltonian. Incidentally, a Hamiltonian of Schneider [32] can be viewed as 
a subsequent singular limit of the Hamiltonian fll.ip . Schneider’s system was mentioned in 
120 ], too, but the relation to van Diejen’s system was not described. 

The original idea behind the present work and [20] was that a natural Poisson-Lie analogue 
of the Hamiltonian reduction treatment [TB] of the BC„ Sutherland system should lead to 
a deformation of this system. It was expected that a special case of van Diejen’s standard 
5-coupling deformation will arise. The expectation has now been conhrmed, although it came 
as a surprise that a singular limit is involved in the connection. 

The outline of the paper is as follows. We start in Section |2] by dehning the reduction 
of interest. In Section E] we observe that several technical results of m can be applied 
for analyzing the reduction at hand, and solve the momentum map constraints by taking 
advantage of this observation. The heart of the paper is Section 01 where we characterize 
the reduced system. In Subsection 14.11 we prove that the reduced phase space is smooth, as 
formulated in Theorem 14.41 Then in Subsection 14.21 we focus on a dense open submanifold 
on which the Hamiltonian fll.ip lives. The demonstration of the Liouville integrability of 
the reduced free flows is given in Subsection 14.31 In particular, we prove the integrability of 
the completion of the system fll.ip carried by the full reduced phase space. Our main result 
is Theorem 14.91 (proved in Subsection 14.4p . which establishes a globally valid model of the 
reduced phase space. We stress that the global structure of the phase space on which the 
flow of fll.ip is complete was not considered previously at all, and will be clarihed as a result 
of our group theoretic interpretation. Section |B] contains our conclusions, further comments 
on the related paper by Marshall I2D1 and a discussion of open problems. The main text is 
complemented by four appendices. Appendix deals with the connection to van Diejen’s 
system; the other 3 appendices contain important details relegated from the main text. 

^We call the limit singular since it involves sending some shifted position variables to infinity. 
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2 Definition of the Hamiltonian reduction 


We below introduce the ‘free’ Hamiltonians and define their reduction. We restrict the 
presentation of this background material to a minimum necessary for understanding our 
work. The conventions follow m, which also contains more details. As a general reference, 
we recommend [7]. 

2.1 The unreduced free Hamiltonians 

We £x a natural numbeiH n > 2 and consider the Lie group SU {2n) equipped with its standard 
quadratic Poisson bracket dehned by the compact form of the Drinfeld-Jimbo classical r- 
matrix, 

^DJ = i Eal3 /\ EjSa, (2.1) 

l<a<0<2n 

where Ea/s is the elementary matrix of size 2n having a single non-zero entry 1 at the af3 
position. In particular, the Poisson brackets of the matrix elements of g E SU( 2 ? 7 ,) obey 
Sklyanin’s formula 

{9 f fi'}su( 2 n) = [9<S)g, roj]. (2.2) 

Thus SU(2n) becomes a Poisson-Lie group, i.e., the multiplication SU(2n) x SU(2n) 
SU(2n) is a Poisson map. The cotangent bundle T*SU(2n) possesses a natural Poisson- 
Lie analogue, the so-called Heisenberg double [33], which is provided by the real Lie group 
SL(2n, C) endowed with a certain symplectic form PP, w. To describe u, we use the Iwasawa 
decomposition and factorize every element K G SL(2n, C) in two alternative ways 

K = = bhOR (2-3) 


with uniquely determined 


gL,gReS\J{2n), 6 ^, 6 ^ G SB(2n). (2.4) 

Here SB( 2 ? 7 ,) stands for the subgroup of SL(2n, C) consisting of upper triangular matrices 
with positive diagonal entries. The symplectic form uj reads 

u = A dgigl^) + ^Atr(d6/j6^^ A dgRg^^). (2.5) 

Before specifying free Hamiltonians on the phase space SL(2n, C), note that any smooth 
function h on SB(2n) corresponds to a function h on the space of positive dehnite Hermitian 
matrices of determinant 1 by the relation 

V 6 G SB(2n). (2.6) 

Then introduce the invariant functions 


C“(SB(2n))®“P”> = {ft e C“(SB(2n)) | /{(ftft*) = { 1 ( 966 * 9 -'), Vg e SU(2ii).6 e SB(2ii)}. 

(2.7) 


These in turn give rise to the following ring of functions on SL( 2 ? 7 ,, C): 


S) = {nE C'~(SL(2n,C)) I n{gLb-^^) = ^( 6 ^), h G C'~(SB(2n))SU(2n)|^ ^2.8) 

•^The n = 1 case would need special treatment and is excluded in order to simplify the presentation. 
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where we utilized the decomposition fl2.3p . An important point is that Sj forms an Abelian 
algebra with respect to the Poisson bracket associated with u fl2.5p . 

The flows of the ‘free’ Hamiltonians contained in Sj can be obtained effortlessly. To 
describe the result, dehne the derivative d^f G C'‘^(SB(2n),5u(2n)) of any real function 
/ G C°°(SB(2n)) by requiring 

f{be^^) = Atr(Xd^/(6)), V6 G SB(2u), VX G Lie(SB(2n)). (2.9) 

s=0 

The Hamiltonian flow generated by "H G through the initial value A'(O) = gL{0)bR{0)~^ is 
in fact given by 

Kit) = 5 (^( 0 ) exp [-td^hibiii0))]b]^^i0), (2.10) 

where 'H and h are related according to fl2.8p . This means that giit) follows the orbit of a 
one-parameter subgroup, while 6_R(t) remains constant. Actually, gnit) also varies along a 
similar orbit, and biit) is constant. 

The constants of motion bi and bn generate a Poisson-Lie symmetry, which allows one to 
dehne Marsden-Weinstein type [I9] reductions. 

2.2 Generalized Marsden-Weinstein reduction 

The free Hamiltonians in Sj are invariant with respect to the action of SU(2n) x SU(2n) on 
SL(2n, C) given by left- and right-multiplications. This is a Poisson-Lie symmetry, which 
means that the corresponding action map 

SU(2n) X SU(2n) x SL(2n, C) ^ SL(2n, C), (2.11) 

operating as 

iVL, Vr, K) ^ VLKg^\ (2.12) 

is a Poisson map. In fl2.11l) the product Poisson structure is taken using the Sklyanin bracket 
on SU(2n) and the Poisson structure on SL(2n, C) associated with the symplectic form u 
02.511 . This Poisson-Lie symmetry admits a momentum map in the sense of Lu [H], given 
explicitly by 

$: SL(2n,C) ^ SB(2n) X SB(2n), <h(A:) = (6^, 6^)- (2.13) 

The key property of the momentum map is represented by the identity 

fie^^Ke-^^) = $>tr(X{/, b^jb^^ + Y{f, bRjb '^), VX, H G 5u(2n), (2.14) 

s=0 

where / G C'°°(SL(2?7,, C)) is an arbitrary real function and the Poisson bracket is the one 
corresponding to u 02 .5p . The map <1> enjoys an equivariance property and one can [TH] 
perform Marsden-Weinstein type reduction in the same way as for usual Hamiltonian actions 
(for which the symmetry group has vanishing Poisson structure). To put it in a nutshell, any 
H E Sj gives rise to a reduced Hamiltonian system by fixing the value of $ and subsequently 
taking quotient with respect to the corresponding isotropy group. The reduced flows can be 
obtained by the standard restriction-projection algorithm, and under favorable circumstances 
the reduced phase space is a smooth symplectic manifold. 

Now, consider the block-diagonal subgroup 

G'+ := S(U(n) X U(n)) < SU(2n). (2.15) 


d 

ds 


A 

ds 
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Since G+ is also a Poisson submanifold of SU(2n), the restriction of fl2.12p yields a Poisson-Lie 
action 

G+xG+x SL(2n, C) ^ SL(2n, C) (2.16) 

of G+ X G+. The momentum map for this action is provided by projecting the original 
momentum map $ as follows. Let us write every element b G SB(2n) in the block-form 


b 


h(l) 6(12) 

On 6(2) 


(2.17) 


and dehne < SB(2n) to be the subgroup for which 6(12) 
denotes the projection 


■6(1) 

6(12)' 

1_V 

•6(1) 

On' 

On 

6(2) _ 

1 —7 

On 

6(2)_ 


On. If tt: SB(2 n) ^ G^ 
(2.18) 


then the momentum map <h+: SL(2?7,, C) —)■ G\ x G'^, is furnished by 


= (ir(()i),!r(()H)). 


(2.19) 


Indeed, it is readily checked that the analogue of fl2.14p holds with X, Y taken from the block- 
diagonal subalgebra of su(2n) and 6^, bji replaced by their projections. The equivariance 
property of this momentum map means that in correspondence to 

K H- with {r]L,VR) e G+ x G+, (2.20) 


one has 


(7r(6i,)7r(6z,)^7r(6i?)7r(6R)^) H- (r7i7r(6z,)7r(6L)'f?7^\r7R7r(6R)7r(6jj)^r7^^). (2.21) 


We briefly mention here that, as the notation suggests, G*^ is itself a Poisson-Lie group 
that can serve as a Poisson dual of G+. The relevant Poisson structure can be obtained by 
identifying the block-diagonal subgroup of SB(2n) with the factor group SB(2n)/L, where L 
is the block-upper-triangular normal subgroup. This factor group inherits a Poisson structure 
from SB(2n), since L is a so-called coisotropic (or ‘admissible’) subgroup of SB(2n) equipped 
with its standard Poisson structure. The projected momentum map is a Poisson map 
with respect to this Poisson structure on the two factors G\ in fl2.19p . The details are not 
indispensable for us. The interested reader may hnd them e.g. in [6]. 

Inspired by the papers [la m EO], we wish to study the particular Marsden-Weinstein 
reduction dehned by imposing the following momentum map constraint: 


^+{K) = /i = (/iL,/iR), where 


e^i/^x) On 
0 e““l 

^ -*-n 


hi? 


e^’ln On 
On e-n 


( 2 . 22 ) 


with some real constants m, v and x. Here, v{x) G SB(n) is the nxn upper triangular matrix 
dehned by 

v{x)jj = 1, v{x)jk = (1 - e"‘^)e^^, j < k, (2.23) 

whose main property is that i/(x)z/(a;)l has the largest possible non-trivial isotropy group 
under conjugation by the elements of SU(n). 

Our principal task is to characterize the reduced phase space 

(2.24) 


7 



















where = {K G SL(2n, C) | = /r} and 


X G 


(2.25) 


is the isotropy group of fi inside G+ x G+. Concretely, is the subgroup of G+ 

consisting of the special unitary matrices of the form 


On ??l(2) 


(2.26) 


where til{2) is arbitrary and 

riL{l)i>{x)v{x)'^riL{l)~^ = z/(a;)z/(a;)^ 


(2.27) 


In words, rjiil) belongs to the little group of i/(a:)z/(a:)l in U(n). We shall see that <I)_,_^(/i) 
and M are smooth manifolds for which the canonical projection 


TT^: <I>+^(/i) M 


(2.28) 


is a smooth submersion. Then M fl2.24p inherits a symplectic form um from u fl2.5p . which 
satishes 

= KM, (2.29) 

where <I)i,i^(/r) —)■ SL(2n, C) denotes the tautological embedding. 


3 Solution of the momentum map constraints 

The description of the reduced phase space requires us to solve the momentum map con¬ 
straints, i.e., we have to find all elements K E <I>+^(/i). Of course, it is enough to do this up 
to the gauge transformations provided by the isotropy group 02.251) . The solution of this 
problem will rely on the auxiliary equation 03.lip below, which is essentially equivalent to the 
momentum map constraint, $+(7^) = /i, and coincides with an equation studied previously 
in great detail in m- Thus we start in the next subsection by deriving this equation. 


3.1 A crucial equation implied by the constraints 

We begin by recalling (e.g. [21]) that any g E SU(2n) can be decomposed as 


9 = 9+ 


cos q 
ising 


ising 
cos g 


h 


+ ) 


where g+, E G+ and g = diag(gi,..., g„) G satishes 

^ > gi > • • • > gn > 0. 


(3.1) 


(3.2) 


The vector g is uniquely determined by g, while g+ and suffer from controlled ambiguities. 

First, apply the above decomposition to gi m K = gLb]^^ £ *h+^(gi) and use the right- 
handed momentum constraint 7r(6ij) = /tr. It is then easily seen that up to gauge transfor¬ 
mations every element of <F^^(/i) can be represented in the following form: 


-1 

O 


cos g i sin g 


e '"In « 

On In 


i sin g cos g 


0 
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(3.3) 



















Here p G SU(n) and a is an n x n complex matrix. By using obvious block-matrix notation, 
we introduce := K 22 and record from fl3.3p that 

= i(sin g)a-t-cos g. (3.4) 


For later purpose we introduce also the polar decomposition of the matrix ff, 


Vt = AT, 


(3.5) 


where T G U(? 7 ,) and the Hermitian, positive semi-dehnite factor A is uniquely determined 
by the relation = A^. 

Second, by writing K = the left-handed momentum constraint vr( 62 .) = Pl tells us 

that bi has the block-form 


bL = 


e“z/(x) 

0 „ 


A 


(3.6) 


with an n X n matrix y. Now we inspect the components of the 2x2 block-matrix identity 


KK* = bi,bi, 


(3.7) 


which results by substituting K from fl3.3p . We hnd that the (22) component of this identity 
is equivalent to 

— e~‘^'"{sm.qY. (3.8) 

On account of the condition fll. 2 p . this uniquely determines A in terms of g, and shows also 
that A is invertible. A further important consequence is that we must have 


qn > 0, (3.9) 

and therefore sing is an invertible diagonal matrix. Indeed, if g„ = 0, then from fl3.4p and 
03.811 we would get (00^)„„ = which is excluded by 01.211 . 

Next, one can check that in the presence of the relations already established, the (12) and 
the (21) components of the identity 03.711 are equivalent to the equation 

y = p(ising)“^[e““ cosg — . (3.10) 

Observe that K uniquely determines g, T and p, and conversely K is uniquely dehned by the 
above relations once g, T and p are found. 

Now one can straightforwardly check by using the above relations that the (11) component 
of the identity 03.7p translates into the following equation: 

p(sing)“^T’^(sing)^T(sing)“^p'^ = v{x)v{x)^. (3-11) 

This is to be satished by g subject to 03.211 . 03.911 and T G U(n), p G SU(n). What makes 
our job relatively easy is that this is the same as equation (5.7) in the paper [TT] by Klimcik 
and one of us. In fact, this equation was analyzed in detail in since it played a crucial 
role in that work, too. The correspondence with the symbols used in m is 

{p,T,smq) {kL,k\^,e^). (3.12) 

This motivates to introduce the variable p G M" in our case, by setting 

singfc = k = l,...,n. (3.13) 
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Notice from fl3.2p and fl3.9p that we have 


0 > Pi > • • • > Pn > —oo. (3-14) 

If the components of p are all different, then we can directly rely on m to establish both the 
allowed range of p and the explicit form of p and T. The statement that pj 7 ^ pk holds for 
j k can be proved by adopting argnments given in |Tll | 12 ] . This proof reqnires combining 
techniqnes of m and na, whose extraction from mm is rather involved. We present 
it in Appendix |Bl otherwise in the next snbsection we proceed by simply stating relevant 
applications of resnlts from mi- 

Remark 3.1. In the context of [m the components of p are not restricted to the half-line 
and both ki and kji vary in U(n). These slight differences do not pose any obstacle to nsing 
the results and techniques of mm- We note that essentially the same equation fl3.1ip 
surfaced in [ 20 ] as well, but the author of that paper refrained from taking advantage of the 
previous analyses of this equation. In fact, some statements of pm are not fully correct. This 
will be specihed (and corrected) in Section |5l 

3.2 Consequences of equation (13.111) 

We start by pointing out the foundation of the whole analysis. For this, we hrst display the 
identity 

iy{x)iy{x)^ = e~^ln + sgn(a;)'0'0^, (3.15) 

which holds with a certain n-component vector v = v{x). By introducing 

w = (3.16) 

and setting p = diag(pi,... ,Pn), we rewrite equation fl3.1ip as 

^ sgn{x)ePww^eP = T-^e^^T. (3.17) 

The equality of the characteristic polynomials of the matrices on the two sides of fl3.17p gives 
a polynomial equation that contains p, the absolute values and a complex indeterminate. 
Utilizing the requirement that \wj\'^ > 0 must hold, one obtains the following result. 

Proposition 3.2. If K given by fl3.3p belongs to the constraint surface <I)+^(p), then the 
vector p fl3.13p is contained in the closed polyhedron 

C := {p e M” I 0 > pi, Pk - Pk+i > \x\/2 {k = l,...,n-l)}. (3.18) 

Proposition 13.21 can be proved by merging the proofs of Lemma 5.2 of [TT] and Theorem 2 of 
pa. This is presented in Appendix |Bl 

The above-mentioned polynomial equality permits to hnd the possible vectors w 03.161) as 
well. If p and w are given, then T is determined by equation 03.17P up to left-multiplication 
by a diagonal matrix and p is determined by fl3.16p up to left-multiplication by elements from 
the little group of v{x). Following this line of reasoning and controlling the ambiguities in 
the same way as in m, one can hnd the explicit form of the most general p and T at any 
hxed p &Cx- In particular, it turns out that the range of the vector p equals Cx- 
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Before presenting the result, we need to prepare some notations. First of all, we pick an 
arbitrary p G Ca, and define the n x n matrix 0{x,p) as follows: 


9{x,p)jk 


sinh (f) 

n 

rr 

■sinh(pj -Pm-fj sinh(pfc - p™ + f )' 

sinh(pfc - pj) 

11 

m=l 

sinh(pj - Pm) sinh(pfc - Pm) 




j^k, (3.19) 


and 


9{x,p)jj:= JJ 


m=l 


sinh(pj -Pm-fj sinh(pj - Pm + f) 
sinh^(pj - Pm) 


(3.20) 


All expressions under square root are non-negative and non-negative square roots are taken. 
Note that 9{x,p) is a real orthogonal matrix of determinant 1 for which 9{x,p)~^ = 9{—x,p) 
holds, too. 

Next, dehne the real vector r{x,p) G M"" with non-negative components 


r{x,p)j = 


1 — e~ 


n 




1 — J-j. y X — 

k=l 

(Mi) 


, j = l,...,n, 


and the real n x n matrix C{x,p), 


C{x,p)aa = r{x,p)a, C{x,p)ij = 5i. 


r{x,p)ir{x,p)j 


l + r{x,p)a 
C,{x,p)ia = -C,{x,p)ai = r{x,p)i, ij 7 ^ a, 

where a = n if x > 0 and a = 1 if x < 0. Introduce also the vector v = v{x): 


(3.21) 


(3.22) 


ln(e^ - 1) _i£ 

v(x]i = \i- -e 2 , j = l,...,n. 


1 — e" 


which is related to v in (13.151) by 


(3.23) 


qUX _ X 

v{x) = \l sgn{x)e~^ —— v{x). 


n 


Finally, define the n x n matrix k(x) as 


r;(x)„ v{x)iv{x)j 

K[x)aa = l^[x)ij = dij -^ ^ , 

y/n n + ^ynv{x)a 

r \ r \ • ■ -L 

k^yXjia — kvyXjai — ,— , 3 / 

In 


(3.24) 


(3.25) 


where, again, a = n if x > 0 and a = 1 if x < 0. It can be shown that both k(x) and C{x,p) 
are orthogonal matrices of determinant 1 for any p ^ Cx- 

Now we can state the main result of this section, whose proof is omitted since it is a direct 
application of the analysis of the solutions of (13.lip presented in Section 5 of [TT] . 
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Proposition 3.3. Take any p G and any diagonal unitary matrix G T„. By using the 
preceding notations define K G SL(2n, C) fl3.3p by setting 

T = e^W{-x,p), p = k{x)C{x,p)-^, (3.26) 


and also applying the eguations fl3.4p . fl3.5p . fl3.8p and fl3.13p . Then the element K belongs 
to the constraint surface and every orbit of the gauge group fl2.25p in 

intersects the set of elements K just constructed. 


Remark 3.4. It is worth spelling out the expression of the element K given by Proposition 
13.31 Indeed, we have 




-1 

O 

s 


Vln - e2P 



e a 

On In 


ie^ 

Vln - 


On e'’ln_ 


using the above dehnitions and 


a = —i 


e'^\/ 6{—x,p) — e"^/e~‘^P — 


(3.27) 


(3.28) 


Remark 3.5. Let us call S the set of the elements K{p, e^^) constructed above, and observe 
that this set is homeomorphic to 


4xT„ = {(p,e''^)} (3.29) 

by its very dehnition. This is not a smooth manifold, because of the presence of the boundary 
of Cx- However, this does not indicate any ‘trouble’ since it is not true (at the boundary of 
Cx) that S intersects every gauge orbit in <I)],]^(/i) in a single point. Indeed, it is instructive 
to verify that if p is the special vertex of Cx for which p^ = {1 — k)\x\/2 for k = 1 ,... ,n, 
then all points K{p,e^^) lie on a single gauge orbit. This, and further inspection, can lead 
to the idea that the variables qj should be identihed with arguments of complex numbers, 
which lose their meaning at the origin that should correspond to the boundary of Cx- Our 
Theorem 14.91 will show that this idea is correct. It is proper to stress that we arrived at such 
idea under the supporting influence of previous works [2Sl[II]. 


4 Characterization of the reduced system 

The smoothness of the reduced phase space and the completeness of the reduced free flows 
follows immediately if we can show that the gauge group acts in such a way on <I)],;^(/r) 
that the isotropy group of every point is just the hnite center of the symmetry group. In 
Subsection 14.11 we prove that the factor of G^ by the center acts freely on $j]^(/i). Then 
in Subsection 14.21 we explain that Cx x T„ provides a model of a dense open subset of the 
reduced phase space by means of the corresponding subset of <I)+^(p) dehned by Proposition 
13.31 Adopting a key calculation from [2U], it turns out that (p, e’”^) G x T„ are Darboux 
coordinates on this dense open subset. In Subsection 14.31 we demonstrate that the reduction 
of the Abelian Poisson algebra of free Hamiltonians fl2.8p yields an integrable system. Finally, 
in Snbsection l4.41 we present a model of the full reduced phase space, which is our main result. 
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4.1 Smoothness of the reduced phase space 

It is clear that the normal subgroup of the full symmetry group G+ x G+ consisting of 
matrices of the form 


( 77 , 77 ) with 77 = diag( 2 ;l„, 2 ;ln), 2 ;^" = ! (4.1) 

acts trivially on the phase space. This subgroup is contained in fl2.25p . The corresponding 
factor group of is called ‘effective gauge group’ and is denoted by G^. We wish to show 
that G^ acts freely on the constraint surface <l)^^(/i). 

We need the following elementary lemmas. 

Lemma 4.1. Suppose that 


cos q i sin o 

9 + ■ ■ 

1 sill q cos q 

with ( 7 +, g'j^, h'^ G G+ and q = diag(gi, ... ,qn) subject to 

^ > qi > ■ ■ ■ > qn > 0- 

Then there exist diagonal matrices mi, m 2 G T„ having the form 


h^= 


9+ 


cos q 1 sm g 
i sin q cos q 


h\ 


(4.2) 


(4.3) 


mi = diag(a,0, 777,2 = diag(6, ,^), ^ G T„_i, a, 6 G Ti, det( 77717772 ) = 1, (4.4) 

for which 

( 5 -+, h+) = (g+diag(777i, 7772 ), dmg{mf^,mf^)h+). (4.5) 

If 04.31) holds with strict inequality | > gi, then mi = m 2 , i.e., a = b. 

Lemma 4.2. Pick any p E Cx and consider the matrix 9{x,p) given by 03.19p and 03.20^ . 
Then the entries 9n^i{x,p) and 9jj-^i{x,p) are all non-zero if x > 0 and the entries 9i^n{x,p) 
and 9j^ij{x,p) are all non-zero if x < 0. 

For convenience, we present the proof of Lemma ITTl in Appendix O The property recorded 
in Lemma 14.21 is known EHUn], and is easily checked by inspection. 

Proposition 4.3. The effective gauge group G^ acts freely on 

Proof. Since every gauge orbit intersects the set S specihed by Proposition 13.31 it is enough 
to show that if {rji, rj^) G maps K E S 03.27P to itself, then {rjL, rjji) equals some element 
( 77 , 77 ) given in 04.ip . For K of the form 03.3p . we can spell out K' = rjiKpf^ as 


Vl{I)p 

0„ ■ 

COS g i sin g 


On 

e ^In r]R{l)ar]R{2) ^ 

o„ 


i sin g cos g 

On 


_ On e^n 


The equality K' = K implies by the uniqueness of the Iwasawa decomposition and Lemma 
10 that we must have 


77l(2) = 77r(1) = 7772, riR{fT)=mi, 77 ^( 1 )^ = p777i, (4.7) 

with some diagonal unitary matrices having the form 04.4p . By using that PRil) = m 2 and 
r]R{2) = mi, the Iwasawa decomposition oi K' = K in 03.27P also entails the relation 

a = m2amf^. (4.8) 
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Because of 03.281) . the off-diagonal components of the matrix equation 04.8p yield 


9{-x,p)jk= ^ k. (4.9) 

This implies by means of Lemma 14.21 and equation 04.4p that mi = m 2 = zln is a scalar 
matrix. But then Pl{^) = follows from pi{l)p = pmi, and the proof is complete. □ 

Proposition 14.31 and the general results gathered in Appendix iDl imply the following theorem, 
which is one of our main results. 

Theorem 4.4. The constraint surface is an embedded submanifold o/SL(2n, C) and 

the reduced phase space M 02.24|) is a smooth manifold for which the natural projection 
TT^: M is a smooth submersion. 

4.2 Model of a dense open subset of the reduced phase space 

Let us denote by S'" C S' the subset of the elements K given by Proposition 13.31 with p in the 
interior of the polyhedron Cx fl3.18p . Explicitly, we have 

= {K{p, I (p, e C, X TJ, (4.10) 

where K{p,e^^) stands for the expression fl3.27p . Note that S'" is in bijection with Cx x T„. 
The next lemma says that no two different point of S'" are gauge equivalent. 

Lemma 4.5. The intersection of any gauge orbit with S'" consists of at most one point. 

Proof. Suppose that 

K’ := K{p\ = plK{p, (4.11) 

with some {pl,Vr) £ G^. By spelling out the gauge transformation as in fl4.6l) . using the 
shorthand sin q = e^, we observe that p' = p since q in fl3.ip does not change under the action 
of G_|_ X G+. Since now we have qi (which is equivalent to 0 > Pi), the arguments applied 
in the proof of Proposition 14.31 permit to translate the equality fl4.1ip into the relations 

Pl( 2) = Pi?(l) = Pi?(2) = m, pL(l)p = pm, (4.12) 

complemented with the condition 

a(p, e^^') = ma(p, e‘'^)m“\ (4.13) 

which is equivalent to 

e^^'9{—x,p) = me^^6{—x,p)m~^. (4.14) 

We stress that m G T„ and notice from fl3.20p that for p ^ Cx all the diagonal entries 
6{—x,p)jj are non-zero. Therefore we conclude from fl4.14p that e^^ = e*'^. This hnishes the 
proof, but of course we can also conhrm that m = zln, consistently with Proposition 14.31 □ 

Now we introduce the map V: SL(2n, C) —)■ M"’ by 

V ■. K = gLbfl ^ p, (4.15) 

dehned by writing g^ in the form fl3.1l) with sing = e^. The map V gives rise to a map 
V : M —)■ M” verifying 

VMK))=V{K), (4.16) 
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where vr^ is the canonical projection fl2.28p . We notice that, since the ‘eigenvalne parameters’ 
Pj (j = are pairwise different for any K G "P is a smooth map. The 

continnity of V implies that 


:= P-\C^) = TT^iSP C M (4.17) 

is an open snbset. The second eqnality is a direct conseqnence of onr foregoing results about 
S and S°. Note that V~^{Cx) = = M. Since is continuous (actually smooth) and 

any point of S is the limit of a sequence in 5°, M° is dense in the reduced phase space M. 
The dense open subset M° can be parametrized by x according to 

(p, e’^) ^ TiPK{p, ep), (4.18) 

which also allows us to view S'® ~ x T„ as a model of M° C M. In principle, the restriction 
of the reduced symplectic form to M° can now be computed by inserting the explicit formula 
K{p,e^P fl3.27p into the Alekseev-Malkin form fl2.5p . In the analogous reduction of the 
Heisenberg double of SU(n, n), Marshall [20] found a nice way to circumvent such a tedious 
calculation. By taking the same route, we have verihed that p and q are Darboux coordinates 
on M°. 

The outcome of the above considerations is summarized by the next theorem. 

Theorem 4.6. M° defined by equation fl4.17p is a dense open subset of the reduced phase 
space M. Parametrizing M° by x according to fl4.18p . the restriction of reduced sym¬ 
plectic form ojm 112.29p to M° is equal to Cj = YPj=i ^ ^Pj Ih.4p . 

4.3 Liouville integrability of the reduced free Hamiltonians 

The Abelian Poisson algebra Sj fl2.8p consists of (G'+ x G_|_)-invariant function^ generating 
complete flows, given explicitly by fl2.10p . on the unreduced phase space. Thus each element 
of Sj descends to a smooth reduced Hamiltonian on M fl2.24p . and generates a complete 
flow via the reduced symplectic form um- This flow is the projection of the corresponding 
unreduced flow, which preserves the constraint surface $j]^(/i). It also follows from the 
construction that Sj gives rise to an Abelian Poisson algebra, Sjm, on Now the 

question is whether the Hamiltonian vector fields of S^m span an n-dimensional subspace of 
the tangent space at the points of a dense open submanifold of M. If yes, then S^m yields a 
Liouville integrable system, since dim(M) = 2n. 

Before settling the above question, let us focus on the Hamiltonian P E Sj dehned by 

n{K) := itr((/<'W)-') = (4.19) 

Using the formula of K{p, e”^) in Remark 13.41 it is readily verihed that 

'H{K{p,eP) = H{p,q-,x,u,v), V(p, e’*^) G x T„, (4.20) 

with the Hamiltonian H displayed in equation fll.ip . Consequently, H in ma is identified 
as the restriction of the reduction of P 04.191) to the dense open submanifold M° 04.17p of 
the reduced phase space, wherein the flow of every element of S)m is complete. 

^More precisely, ^ = C'°°(SL(2n, C))SU( 2 n)xSU( 2 «)^ 
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Turning to the demonstration of Liouville integrability, consider the n functions 


W/=(A') := = r..., 


n. 


(4.21) 


The restriction of the corresponding elements of S)m on M° ~ x gives the functions 


Hk{p,q) = ^tr 


-e^a 


“I k 


(4.22) 


where a has the form fl3.28p . These are real-analytic functions on x T„. It is enough 
to show that their exterior derivatives are linearly independent on a dense open subset of 
Cj; X T„. This follows if we show that the function 


f{p, q) = det [dqHi, dqH 2 ,d^jHn] (4.23) 

is not identically zero on x T^- Indeed, since / is an analytic function and x T„ is 
connected, if / is not identically zero then its zero set cannot contain any accumulation point. 
This, in turn, implies that / is non-zero on a dense open subset of x T„ ~ M°, which is 
also dense and open in the full reduced phase space M. In other words, the reductions of Hk 
{k = 1,... ,n) possess the property of Liouville integrability. It is rather obvious that the 
function / is not identically zero, since involves dependence on q through and lower 
powers of It is not difficult to inspect the function f{p,q) in the ‘asymptotic domain’ 
where all differences \pj — pm\ (^ 7^ j) tend to infinity, since in this domain a becomes close 
to a diagonal matrix. We omit the details of this inspection, whereby we checked that / is 
indeed not identically zero. 

The above arguments prove the Liouville integrability of the reduced free Hamiltonians, 
i.e., the elements of S)m- Presumably, there exists a dual set of integrable many-body Hamil¬ 
tonians that live on the space of action-angle variables of the Hamiltonians in S)m- The 
construction of such dual Hamiltonians is an interesting task for the future, which will be 
further commented upon in Section [5l 

4.4 The global structure of the reduced phase space 

We here construct a global cross-section of the gauge orbits in the constraint surface (/r). 
This engenders a symplectic diffeomorphism between the reduced phase space {M,um) and 
the manifold {McUc) below. It is worth noting that {McUc) is symplectomorphic to 
carrying the standard Darboux 2-form, and one can easily find an explicit symplectomor- 
phism if desired. Our construction was inspired by the previous papers Eattn, but detailed 
inspection of the specific example was also required for finding the final result given by The¬ 
orem SSI After a cursory glance, the reader is advised to go directly to this theorem and 
follow the definitions backwards as becomes necessary. See also Remark 14.101 for the rationale 
behind the subsequent definitions. 

To begin, consider the product manifold 

■= C”-i X D, (4.24) 

where D stands for the open unit disk, i.e., D := {w G C : |w| < 1}, and equip it with the 
symplectic form 

E J A .J- I "^dzn/\ dzn X 

dzj A dzj -|- —-^—. (4-25) 
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The subscript c refers to ‘complex variables’. Define the surjective map 


Z^ \ Me, {P, e*'') i-t z{p, e’'') 

by the formulae 

1 ^ 

Zj{P,e^^) = {pj-pj+i-\x\/2)2 JJ j = 

k=j+l 

1 

fc=i 


(4.26) 


(4.27) 


Notice that the restriction Zx of Z^ to Cx x T„ is a diffeomorphism onto the dense open 
submanifold 

n 



n 

(4.28) 

It verihes 



Z*(u^) = uj = ^dqj A dpj, 

(4.29) 

which means that Z^ 
^ Cx X 

c is a symplectic embedding of {Cx x T„,d)) into {Mc,ujc). 
operates according to 

The inverse 


Pi{z) = log(l - |x„n, pj{z) = log(l - 



ZnZl 


ZnZll 







- 1^1/2) (J = 2 ,..., n) 

k=l 

(m = 2,..., n — 1), e’'^"(x) = . 

\Zn-l\ 


(4.30) 


It is important to remark that the pk{z) (/c = 1,... ,n) given above yield smooth functions 
on the whole of while the angles are of course not well-dehned on the complementary 
locus of M°. Our construction of the global cross-section will rely on the building blocks 
collected in the following long dehnition. 

Definition 4.7. For any {zi ,..., Zn-i) G consider the smooth functions 


Qjk{x, z) 


■sinh(X)J=i + (^ - i)kl/2 -xjZ)' 

_ sinh(X;f=j zizi + (A; - j)\x\/2) 


1 

2 


1 < j < A; < n, 


(4.31) 


and set Qjk{x, z) 


■= Qkj{—x, z) for j > k. Applying these as well as the real analytic function 


Ay) ■■= 


I sinh(|/) 

y 


y J(0) := 1, 


(4.32) 
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and recalling fl3.2ip . introdnce the n x n matrix ({x, z) by the formnlae 


C{x,z)aa = r{x,p{z))a, , z)aj = -({x, z)ja, j 7 ^ a, 

C(x, z)jn = 


sinh(|) ZjJ(zjZj) 


sinh(^) sinh(zjZj + 


JJ Qj((x,z), x> 0 , j^n, 


2 ^ r=i 


((x,z)ji = 


sinh(|) Zj-iJ(zj-iZj-i) 


sinh(^) sinh(zj-iZj-i — |) 


n x< 0 , j^l, 


(4.33) 


2 ' 1=1 


N _ r , i{^^Aial{x,z)ak . , , 
C{x, Zjjk — Oj^k + ^ . 1/ N ) Ji^ 7^ 


1 + C{x,z)c 

where a = n if x > 0 and a = 1 if x < 0. Then introdnce the matrix 6 {x, z) for x > 0 as 


9{x, z)jk = 


sinh(^)sgn(fc - j - l)C{x, z)jnCi-x, z)ik 
sinh(Er=io!fc)' zeze+\k-j - 1 \^) 
sinh(f) 


, fc 7^ j + 1 , 


= - 


smh.{zjZj + |) 




(4.34) 


t=i 


and for X < 0 as 

9{x, z) = 9{—x, z)^. 

Finally, for any G Me dehne the matrix j{x, z) = diag( 7 i,..., %) with 


(4.35) 


7iz)i = ZnV‘2 - ZnZn, 7(x,z)j = 


1 - (1 - 


, j = 2,...,n, (4.36) 


and the matrix 

d(x, M, n, z) = —i — e~^^lnj 9(—x, z) — e'^e~^^^^j'(x, z)^], (4.37) 

nsing the constants x,u,v snbject to fll. 21 ) . 

Althongh the variable Zn appears only in 71 , we can regard all objects dehned above as 
smooth fnnetions on Me, and we shall do so below. 

The key properties of the matrices (, 9, a and 7 are given by the following lemma, which 
can be verihed by straightforward inspection. The role of these identities and their origin 
will be enlightened by Theorem 14.91 

Lemma 4.8. Prepare the notations 


X(x) '■ 

:= diag(r 2 ,..., 

.Tn, 1) 

if X > 0 

and 

X(x) 

:= diag(l,ri,.. 

. .,r„_i) if X < 0, 

(4.38) 

X{x) 

with 

:= diag(l,r 2 ,. 

• • ) Ai) 

n 

if X > 0 

and 

X{x) 

:= diag(ri,... 

n 

,r„_i,l) if x<0 

(4.39) 


= 

k=j 

if X > 0 

and 


= Y[ 

k=j 

if X < 0. 

(4.40) 
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Then the following identities hold for all {p, e^^) E x T„.' 

C{x,z{p,e'^)) = r(^)C(a;,p)r(“J, (4.41) 

9{x, z{p, e'^)) = T^^)9{x,p)ffJ^, (4.42) 

7 (x, z{p, e'^)) = (4.43) 

a{x, u, V, z{p, e"^)) = e~^'^^^)a{x, u, v, p, e''^)r("J. (4.44) 


Here we use Definition \f.'l\ and the functions on x introduced in Subsection \3.^ 

For the verification of the above identities, we remark that the vector r fl3.2ip can be 
expressed as a smooth function of the complex variables as 


r{x,p{z))j = 


sinh( 


sinh( 


Yl Qjk{x,z), j = l,...,n. 


(4.45) 


2 ' k=l 

{k^j) 


With all necessary preparations now done, we state the main new result of the paper. 
Theorem 4.9. The image of the the smooth map K : Me —^ SL(2n, C) given by the formula 


K(z) = 


1 

s 

k 

o 

1 

o 

3 


^{x, z) iePP 


e "In a{x,u,v, z) 

On In 


_ 1 


0 e"l 

'-'n ^ ■‘■n 


(4.46) 


lies in <F+^(/r), intersects every gauge orbit in precisely one point, and K is injective. The pull¬ 
back of the Alekseev-Malkin 2-form u fl2.5p by K is Uc 114.25p . Conseguently, ti^o K: Me —)■ 
M is a symplectomorphism, whereby {Me, (Xe) provides a model of the reduced phase space 
{M,ojm) defined in Subsection \2.2\ 


Proof. The proof is based upon the identity 


K{z{p,P'^)) = 


K{x)T(^,e:)K{x) 


-1 


0 . 


0 . 

pxy 


3-1? 


K{p,e 




r(x)e 0 ^ 

Pp. 


-\ -1 


0 . 


, V(p, 6 “^) eCxX T„, 
(4.47) 

which is readily seen to be equivalent to the set of identities displayed in Lemma ITBl It means 
that K{z{p, e''^)) is a gauge transform of K{p, e'^) in fl3.27p . Indeed, the above transformation 
of K{p, e''^) has the form fl 2 . 20 p with 


Vl = c 


lz{x)Ti^x)t^{x) 


-1 


0 . 


0 . 

px)G 


-iq 


Vr = c 


r(x)e On 

On '^(.D 


(4.48) 


where c is a harmless scalar inserted to ensure det{pL) = det{pR) = 1. Using fl3.25p and 
fl4.38p . one can check that k{x)t(^x)R{x)~^v{x) = v{x) for the vector v{x) in fl3.24p . which 
implies via the relation fl3.15p that {riL,riR) belongs to the isotropy group 112.25p . the 
gauge group acting on $],]^(p). 

It follows from Proposition 13.31 and the identity 114. 4 7p that the set 

S := {k{z) \ zEMe} (4.49) 

lies in $^^(/i) and intersects every gauge orbit. Since the dense subset 


:= {k{z) I z e M°} 


(4.50) 
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is gauge equivalent to S° in fl4.10p . we obtain the equality 


k*{u)=u; (4.51) 

by using Theorem 14.61 and equation fl4.29p . More precisely, we here also utilized that K*{u) 
is (obviously) smooth and M° is dense in M^. 

The only statements that remain to be proved are that the intersection of S with any 
gauge orbit consists of a single point and that K is injective. (These are already clear for 
S° <Z S and for Now suppose that 


kiz') 


VLik 


Or 


Vl{ 2 ) 


K{z) 


vkk On 

On Vr{2) 


(4.52) 


for some gauge transformation and z, z' G Mf~. Let us observe from the dehnitions that we 
can write 

7 ( 0 :, z) 

igPh) zY 

where sing( 2 ;) = with | > gi > • • • > gn > 0, and D{z) is a diagonal unitary matrix 
of the form D{z) = diag((ii, ln-i,di, In-i)- Then the uniqueness properties of the Iwasawa 
decomposition of SL(2n, C) and the generalized Cartan decomposition fl3.ip of SU(2n) allow 
to establish the following consequences of fl4.52p . First, 


= D(z) 


COS q{z) ising( 2 ;) 
ising(z) cosg( 2 :) 


D(z), 


(4.53) 


p(z) = p{z'). 


Second, using Lemma ITTl 



0 „ ■ 


m 2 On 

_ 0 , 

hi?(2). 


On rni_ 


for some diagonal unitary matrices of the form 04.41) . Third, we have 

aiz') = r]R{l)a{z)T]R{2)~^ = m2a{z)mk ■ 


(4.54) 

(4.55) 


(4.56) 


For dehniteness, let us focus on the case a; > 0. Then we see from the dehnitions that the 
components dfc+yfc and ai^n depend only on p(z) and are non-zero. By using this, we hnd 
from 04.56P that mi = m 2 = Cln with a scalar C, and therefore 

6 i{z') = a{z). (4-57) 

Inspection of the components (1, 2),..., (1, n — 1) of this matrix equality and 04.54p permit 
to conclude that Z 2 = = Zn-i, respectively. Then, the equality of the (2,?7,) 

entries in 04.57p gives z'^ = zi which used in the (1,1) position implies z'^^ = Zn- Thus we see 
that z' = z and the proof is complete. (Everything written below fl4.56p is quite similar for 

X < 0.) □ 

Remark 4.10. Let us hint at the way the global structure was found. The hrst point to 
notice was that all or some of the phases e'h cannot encode gauge invariant quantities if 
p belongs to the boundary of Cx, as was already mentioned in Remark 13.51 Motivated by 
HU, then we searched for complex variables by requiring that a suitable gauge transform of 
K[p, in 03.271) should be expressible as a smooth function of those variables. Given the 
similarities to HU, only the dehnition of Zn was a true open question. After trial and error, 
the idea came in a flash that the gauge transformation at issue should be constructed from 
a transformation that appears in Lemma 1C.11 Then it was not difficult to hnd the correct 
result. 
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Remark 4.11. Let us elaborate on how the trajectories p{t) corresponding to the flows of 
the reduced free Hamiltonians, arising from Tik fl4.2ip for k = can be obtained. 

Recall that for k = 1 the reduction of "Hi completes the main Hamiltonian H fll.ip . Since 
HkiK) = hkibfi) with hk{b) = ^tr(6l6)*', the free flow generated by Tik through the initial 
value R'(O) = fi'L(0)6jj^(0) is given by 02.1011 with d^hk{b) = i{b^bY. Thus the curve giit) 
O2.10p has the form 

= 5 'l(0 ) exp(-h£(0)*') with £(0) = 6ij(0)^6R(0). (4.58) 


The reduced flow results by the usual projection algorithm. This starts by picking an initial 
value 2:(0) G Me and setting i^(0) = K{z{0)) by applying 04.4611 . which directly determines 
(7 l(0 ) and 5i?(0) as well. Then the map V 04.15p gives rise to p{t) via the decomposition of 
giit) G SU(2n) as displayed in 03.111 . that is 


p(t) = nm)- 


(4.59) 


More explicitly, if Vit) stands for the (11) block of giii)) then the eigenvalues of V{t)V{ty 
are 

(T{V{t)V{t)^) = {cos^ qj{t) I j = 1,..., n}, (4.60) 

from which Pj{t) can be obtained using 03.1311 . In particular, the ‘particle positions’ evolve 
according to an ‘eigenvalue dynamics’ similarly to other many-body systems. This involves 
the one-parameter group , where /I(0) is the initial value of the Lax matrix (cf. 04.22p j 


—e’’Q;(^) 

—e'"a{zy (e“^'"l„ -h a{zya{z)) ’ 


(4.61) 


where we suppressed the dependence of a 04.37P on the parameters x,u,v. A more detailed 
characterization of the dynamics will be provided elsewhere. 


5 Discussion and outlook on open problems 

In this paper we derived a deformation of the trigonometric BC„ Sutherland system by means 
of Hamiltonian reduction of a free system on the Heisenberg double of SU(2n). Our main 
result is the global characterization of the reduced phase space given by Theorem 14.91 The 
Liouville integrability of our system holds on this phase space, wherein the reduced free flows 
are complete. These flows can be obtained by the usual projection method applied to the 
original free flows described in Section [2j 

The local form of our reduced ‘main Hamiltonian ’ dni) is similar to the Hamiltonian 
derived in [20], which deforms the hyperbolic BC„ Sutherland system. However, besides a 
sign difference corresponding to the difference of the undeformed Hamiltonians, the local 
domain of our system, Cx x T„ in fll.3p . is different from the local domain appearing in [20], 
which in effect has the form Cj, x with the open polyhedron^ 


Cj, := {p G M*" I Pn > 0, Pk - Pk+i > \x\/2 {k = l,...,n-l)}. (5.1) 

We here wish to point out that Cj, x T„ is not the full reduced phase space that arises from 
the reduction considered in [20]. In fact, similarly to our case, the constraint surface contains 
a submanifold of the form Cj, x T„ in the case of [20], where Cj, is the closure of Cj,. Then a 

^The notational correspondence between [20] and the present paper is: a, x, y) O (p, q, , e~'", e““). 
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global model of the reduced phase space can be constructed by introducing complex variables 
suitably accommodating the procedure that we utilized in Subsection 14.41 Erroneously, in 
123 the full phase space was claimed to be x T„; the details of the correct description will 
be presented elsewhere. 

Throughout the text we assumed that n > 1, but we now note that the reduced system 
can be specialized to n = 1 and the reduction procedure works in this case as well. The 
assumption was made merely to save words. The formalism actually simplihes for n = 1 
since the Poisson structure on (7+ = S(U(1) x U(l)) < SU(2) is trivial. 

As explained in Appendix I3 the Hamiltonian fll.ljl is a singular limit of a specialization 
of the trigonometric van Diejen Hamiltonian |35], which (in addition to the deformation 
parameter) contains 5 coupling constants. As a result, at least classically, van Diejen’s 
system can be degenerated into the trigonometric BC„ Sutherland system either directly, as 
described in [35], or in a roundabout way, going through our system. Of course, a similar 
statement holds in relation to hyperbolic BCn Sutherland and the system of 123. 

Except in the rational limit [26], no Lax matrix is known that would generate van Diejen’s 
commuting Hamiltonians. In the reduction approach a Lax matrix arises automatically, in 
our case it features in equations fl4.22p and fl4.6ip . This might be helpful in searching for 
a Lax matrix behind van Diejen’s 5-coupling Hamiltonian. The search would be easy if 
one could derive van Diejen’s system by Hamiltonian reduction. It is a long standing open 
problem to find such derivation. Perhaps one should consider some ‘classical analogue’ of the 
quantum group interpretation of the Koornwinder (BC„ Macdonald) polynomials found in 
[25] . since those polynomials diagonalize van Diejen’s quantized Hamiltonians [27] . 

Another open problem is to construct action-angle duals of the deformed BC„ Sutherland 
systems. Duality relations are not only intriguing on their own right, but are also very useful 
for extracting information about the dynamics [28l |29l [301 EZ]. The duality was used in 
mm to show that the hyperbolic BC„ Sutherland system is maximally superintegrable, the 
trigonometric BC„ Sutherland system has precisely n constants of motion, and the relevant 
dual systems are maximally superintegrable in both cases. These studies, which were heavily 
influenced by Pusztai’s paper m (see also ini). may provide inspiration for a future inves¬ 
tigation of the dualities for the deformed BC„ Sutherland systems. We here only remark 
that deformed dual systems should arise from considering the reduction of alternative sets of 
commuting free Hamiltonians on the pertinent Heisenberg doubles. 

After we finished our work, there appeared a preprint [22] dealing with the quantum 
mechanics of a lattice version of a 4-parameter Inozemtsev type limit of van Diejen’s trigono¬ 
metric/hyperbolic system. The systems studied in [20] and in our paper correspond to further 
limits of specializations of this one. The statements about quantum mechanical dualities con¬ 
tained in [38] and its references should be related to classical dualities. 

We hope be able to return to some of these questions in the future. 

Acknowledgements. L.F. wishes to thank S. Ruijsenaars for suggesting the role of the 
kind of singular limit described in Appendix [^ This work was supported in part by the 
Hungarian Scientihc Research Fund (OTKA) under the grant K-111697. The work was 
also partially supported by the European Union and the European Social Fund through the 
project TAMOP-4.2.2.D-15/1 at the University of Szeged. 
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A Links to systems of van Diejen and Schneider 

Recall that the trigonometric BC„ van Diejen system [35] has the Hamiltonian 

n 

//.D(A,e) = ^ (cosh(«,)Vj(A)‘''^V_,(A)‘/" - |Vj(A) + V^,(A)]/2), 

i=i 

with \/±j {j = 1 ,... ,n) defined by 

n 

V±j(A) = w(±Aj) v(±Aj + Afc)v(±Aj - A^), 
k=l 

and V, w denoting the trigonometric potentials 

, , sin(/i + z) , , , sin(/ro + z) cosifii + z) sinf/i'o + z) cosf/rj + z) 

V ^ = — ■ / X and w ^ =-—--—-—---—, 

sm( 2 ;j sm(zj cos( 2 ;j sm( 2 ;j cos( 2 ;j 

where fi, fii, are arbitrary parameters. By making the substitutions 

Aj—ti(pj + R), . , po —)■ i(5'o +-R), t if?! + 7r/2, 

a Vj and /i-> i^/2 , da ' ^ , /o 

/^o H5'o ~ /^i ifi'i + 

the potentials become hyperbolic functions and their R ^ oo limit exists, namely 
lim v(±(Aj + Afc)) = lim v(±(Aj - Afc)) = Pk) y ■ ^ 

R^-oo ^ ^ ^ ^ smh{pj-pk) 


(A.l) 


(A. 2 ) 


(A.3) 


(A.4) 


and 


lim w(±Aj) = e®'" — e" 

R—^oo 


= p9o-go±{9l+9i)-^Pj _ pM90+9o+9l+9'i) ^ \/j_ 

In the 1 -particle case we have V±(A) = w(±A), thus H^b takes the following form 
hfvD(A,6') = cosh(6')w(A)^Aw(—— [w(A) -|-w(—A)]/ 2 . 

By utilizing flA.bp one obtains 

lim w(A)1Aw(-A)B 2 = h - {e^ 9 o + e-^g'oy-^P + ^250-25^#! 

C?_h -1 


(A.5) 

(A. 6 ) 

(A.7) 


R—^oo 

lim [w(A) -|-w(—A)]/2 = 

R -^00 


^9o-9o+9l+9{ _|_ q9o-9o-9i-9'i 


- 2 p 


cosh{go + Pq + gi + g[). 


(A. 8 ) 


Equating the i? —>■ 00 limit of i7vD(A,6*) (I A. 71) with the Hamiltonian H{p,q;x,u,v) fll.ip 
yields a system of linear equations involving go, gi, g^, g[ as unknowns and u, v as parameters. 
Actually, four sets of linear equations can be constructed, each with infinitely many solutions 
depending on one (real) parameter, but these sets are ‘equivalent’ under the exchanges: 
9 o 9 o or 9 i ^ 9 i- Therefore it is sufficient to give only one set of solutions, e.g. 

go = v-u, 5'o = 0, gi = u + v-g[, g[eR. (A.9) 

Setting g = X and g[ = 0 provides the following special choice of couplings in flA.4p 

p = ix/ 2 , pQ = i{v — u + R), p'q =—iR, pi = i{u + v) + 7 i/ 2 , p[ = 7 i/ 2 , (A.10) 
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and one finds the following 

lim //vd(A(p,-R), 6 *(g)) =—H{p,q;x,u,v) + cosh (2u). 


R^oo 


(A.ll) 


In the n-particle case, by nsing flA.Sp and flA. 6 D it can be shown that with flA.lOp one has 


lim H^D[X{p,R),9{q)) = -H{p,q;x,u,v) + cosh ((j - l)x + 2 n), (A. 12 ) 

i=i 


i.e., the Hamiltonian H fll.ll) is recovered as a singnlar limit of i^vD HA.ip . 

Consider now the fnnction H{p, g; x, u, v) and introdnce the real parameter a throngh the 
snbstitntions 

u^u — a, v^v — a (A.13) 

and apply the canonical transformation 

Pj —Qj + cr, Qj ^ —-Pj, Vj. (A. 14) 


Then we have 

lim H{p{Q, a), g(P), x, u{a),v{a)) = ifsch(Q, P, x, n), 

cr—yoo 

with Schneider’s [32] Hamiltonian 


Hschi^Qy P) X, n) 



j=i 


^cos(P,) Yl 

j=l k=l 


sinh^ (I) ■ 

sinh^(Qj -Qk). 


1 

2 


(A.15) 


(A.16) 


Remark A.l. (z) In flA.4D only two of the fonr external held conplings po, /Zg, pi, p[ are scaled 
with R. However, scaling all fonr of these parameters also leads to an integrable Rnijsenaars- 
Schneider type system with a more general 4-parameter external held. For details, see Section 
H.B of [36] • {ii) The connection to Schneider’s Hamiltonian was mentioned in Remark 7.1 
of [20] as well, where a singnlar limit, similar to flA.lSp was taken. 


B Proof of a key result 


In this appendix we prove Proposition l3. 21 which states that the range of the ‘position variable’ 
p is contained in the closed thick-walled Weyl chamber fl3.18p . 


Proof of Proposition \3. gt According to fl3.17p the matrices and + sg'a.{x)e^ww^eP 

are similar and therefore have the same characteristic polynomial. This gives the identity 


- A) = 

i=i i=i 


A) -h sgn(x) ^ 



n 


2 Yl 

k=l 

ik¥=j) 



(B.l) 


where A is an arbitrary complex parameter. The constraint on p arises from the fact that 
I'iXmP {m = 1,... ,n) mnst be non-negative and not all zero becanse of the dehnition fl3.16p . 

Let ns assnme for a moment that the components of p are distinct snch that pi > ■ ■ ■ > pn- 
This enables ns to express for all m G {1,..., n} from the above eqnation by evalnating 
it at n different valnes of A, viz. A = m = 1 ,..., n. We obtain the following 


g2pj-|-X _ ^2pm 

n ;™|2 = sgn(x)(l-e-") J] m = l,...,n. 

i=i 


(B.2) 
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For a: > 0 and any p with pi > ■ ■ ■ > Pn the formnla (IB.2p implies that > 0 and for 

m = 1,... ,n — 1 we have \wm\^ > 0 if and only if pm — Pm+i > x/2. Similarly, if x < 0 
and p G M” with pi > ■ ■ ■ > Pm then flB.2l) implies > 0 and for m = 2,... ,n we have 
> 0 if and only if Pm-i — Pm > —x/2. In snmmary, if pi > • • • > pn, then \wm\'^ > 0 
Vm implies that p & Cx. 

Now, let ns prove onr assnmption, that all components of p must be different. Indirectly, 
suppose that some (or maybe all) of the pjS coincide. This can be captured by a partition 
of the positive integer 

n = ki + ■ ■ ■ + kr, (B.3) 

where r < n (or equivalently, at least one integer ki,... ,kr must be greater than 1) and the 
indirect assumption can be written as 


Pi ■ ■ ■ Pkx) Pki+1 ■ ■ ■ Pk\+k2-i ■ ■ ■ I Pk\-\ \-k, — i + l ' ' ' Pk\-\ - \-kr — Pn- (B.d) 

Then fIB.ip can be reformulated as 


n(A, - - Xp + sgn(x) Z^A^A^e-^ - X)^-'^ J] (A,e- 

j=l j=l m=l j=l 


where we introduced r distinct variables 

Ai = A2 = ..., Ar = 

and r non-negative real variables 

Zi = \wp -|- • • • -f \wkp^ Z 2 = Iwfci + lp + • • • + |Wfci+fc2p, 
. . . , Zj. -—1 + 11 T ■ ■ ■ T . 



(B.5) 

(B.6) 


(B.7) 


Notice that Zi + ■ ■ ■ + Z^ = = sgn(x)e ^(e"'^ — 1) > 0, therefore at least one of the Zp 

must be positive. Next, we dehne the rational function of A 


Q(A,x, A) 


n (A, - A)^^' 

{AjC-^ - A)^7-i’ 


(B. 8 ) 


and use it to rewrite fIB.Sp as 

r r r 

Q(A,x, A) = JJ(Aje""' - A) -h sgn(x) ^ Z^A^ JJ (A^e"^ - A). (B.9) 

j=l m=l j=l 

{j¥=m) 

The above equation implies that all poles of Q are apparent, i.e., there must be cancelling 
factors in its numerator. This observation has a straightforward implication on the A’s. 

(+) For every index m 6 {1,..., r} with km > 1, there exists an index s G {1,..., r} s.t. 

As = Am^~^ and ks> km — 

The quantities Zm = Zm{A,x) can be uniquely determined by evaluating flB.Qp at r different 
values of the parameter A, namely A^ = A^e”^ (m = 1,..., r). However, there are 3 disjoint 
cases which are to be handled separately. 
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Case 1: = 1 and G {1,..., r}: = A^e Then we find 


= sgn(a;)(l - " ) ' > 0- (B.IO) 

j=l \ 3 m / 

Case 2: km > ^ and kg = km — 1- Then we find 

Z™ = (-l)'=-+'sgn(a:)(l-e-")e("-^™)" f[ > 0. (B.ll) 

j=l \ 3 m J 

(ji^rn.s) 


Case 3: = 1 and 3s G {1,..., r}: A^ = A^e * or /c^ > 1 and kg > km — 1. Then we get 


Zm — 0. 


(B.12) 


Since there is at least one Zm which is positive, the set of indices belonging to Case 1 or 
Case 2 must be non-empty. Introduce a real positive parameter e and associate to every 
degenerate configuration flB.4H a continuous family of configurations, denoted by p{e), with 
components p{s)i,... ,p{e)n defined by the formulae 


exp(2p(£)a + ae) = Ai, a = l,...,fci, 
exp(2p(£) = a = 1, • • •, kj, j = 2,..., r. 

This way coinciding components of p flB.4p are ‘pulled apart’ to points successively separated 
by e/2. It is clear that with sufficiently small separation the configuration p{e) sits in the 
chamber {f G M"' | 0 > fi > • • • > f„}. For such non-degenerate configurations p{e), let us 
consider the expressions 


\we{p{e),x)\‘^ 


sgn(a;)(l 


-on 

3 = 3 

im) 


^2p(e)j+x _ ^2p{e)e 
Q‘3p(e)j _ Q2p{e)e ’ 


= 1 , 


,n. 


(B.14) 


which give the unique solution of equation fIB.ip at p{e). The limits \\m.e^o\w(,{p{e),x)\^ 
exist, and do not vanish for £ = fci if belongs to Case 1 or Case 2. For such 

I = ki + ■■■ + km must have 


\im\wk,+...+k^{p{e), x)\^ = Z^(A,a;) > 0, (B.15) 

£->■0 


where Zm is given by fIB.lOp in Case 1 and by fIB.llI) in Case 2. It can be also seen that 


{ i ^ {fcl, ki ^2, . . . , ki kr} 

or (B.16) 

i = ki + ■■■ + km with km from Case 3, 

i.e., \we{p{e),x)\‘^ vanishes identically except for the components in flB.lSp . Notice that 
for a small enough e some coordinates of p{e) are separated by less than \x\/2. Thus, as 
it was shown at beginning the proof, we have |tc^(p(£),a;)p < 0 for some index £, which 
might depend on e. Moreover, flB.1611 implies that the index in question must have the form 
^ = ki + ■ ■ ■ + km* for some m* appearing in flB.15F But since the number of indices is 


26 

























finite, a monotonically decreasing sequence tending to zero can be chosen such that 

|tCfci+---+fc^* (p(£Ar), a:)p < 0 for all N. This together with flB.16D gives the contradiction 

0> lim |wfc^+...+fc (p(£jv),a;)p = Z^*(A,a;) > 0 (B.17) 

A'—>-oo 

proving that all components of p must be distinct. This concludes the proof. □ 

The above proof is a straightforward adaptation of the proofs of Lemma 5.2 of m and 
Theorem 2 of na. We presented it since it could be awkward to extract the arguments from 
those lengthy papers, and also our notations and the ranges of our variables are different. 


C Proof of an elementary lemma 

We here prove the following equivalent formulation of Lemma 14.11 
Lemma C.l. Suppose that ^ > qi > ■ ■ ■ > Qn > 0 and 



0„ ■ 

COS q i sin q 


On 


cos q i sin g 

_ On 


i sin q cos q 




i sin g cos g 


for pl^Vr ^ G+- Then 

hn(l) = hR(2) = mi, = pr{ 1) = m2 (C.2) 

with some diagonal matrices mi, m 2 G having the form 

mi = diag(a,^), m 2 = diag(6,^), ^ G T„_i, a,6GTi, det(mim 2 ) = 1. (C.3) 


If in addition \> qi, then mi = m 2 - 

Proof. The block off-diagonal components of the equality flC.ip give 

77 l( 1 ) = (sing) 77 ij( 2 )(sing)“\ 77 ^( 2 ) = (sin g) 77 u(l) (sin (C.4) 

Since , the hrst of these relations implies Pr{ 2) = {smq)^pR{2){smq)~^. As 

the entries of (sing) are all different, this entails that Pr{2) is diagonal, and consequently we 
obtain the relations in flC.2p with some diagonal matrices mi and m 2 . On the other hand, 
the block-diagonal components of flC.ip require that 

cosq = pLiI){cosq)pR{l)~^, cosq = pL{2){cosq)pR{2)~^. (C.5) 


Since cosg^ 7 ^ 0 for fc = 2,..., n, the formula fIC.Sp follows. If an addition cosgi 7 ^ 0, then 
we also obtain from fIC.Sp that a = b, i.e., mi = m 2 = m with some m G T„. □ 


D Auxiliary material on Poisson-Lie symmetry 


The statements presented here are direct analogues of well-known results mm about Hamil¬ 
tonian group actions with zero Poisson bracket on the symmetry group. They are surely 
familiar to experts, although we could not hnd them in a reference. 

Let us consider a Poisson-Lie group G with dual group G* and a symplectic manifold P 
equipped with a left Poisson action of G. Essentially following Lu [18] (cf. Remark [D.4p . we 
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say that the G-action admits the momentum map P —)■ G* if for any X E Q, the Lie 
algebra of G, and any / E C°°{P) we have 


m = (X,{/,V^}(p)^(p)-i), Vp G P, (D.l) 


where Xp is the vector held on P corresponding to X, stands for the canonical pairing 
between the Lie algebras of G and G*, and the notation pretends that G* is a matrix group. 
Using the Hamiltonian vector held Vf dehned by Cy^h = —{/, h} (Vh E C°°{P)), we can 
spell out equation fID.ip equivalently as 

(£x,/)(p) = -(V(C*(rtfi*(rt-.)((r>pV')(V/(p)))), VpeP, (D.2) 


where Dp^|J: TpP -E- T^^p^G* is the derivative, and denotes the right-translation on 

G* by Since the vectors of the form Vf{p) span TpP, we obtain the following charac¬ 

terization of the Lie algebra of the isotropy subgroup Gp < G of p G P. 

Lemma D.l. With the above notations, we have 

Lie(Gp) = [(P,/,(p)P^(p)-i)(lm(PpV’))]'^. (D.3) 


This directly leads to the next statement. 

Corollary D.2. An element fi E G* is a regular value of the momentum map if and only 
z/Lie(Gp) = {0} for every p E = {p G P | '0(p) = h}- 

Let us further suppose that ip: P —)■ G* is G-equivariant, with respect to the appropriate 
dressing action of G on G*. Then we have 

Gp<Gp, ypEij-\p). (D.4) 


Here Gp and Gp refer to the respective actions of G on P and on G*. Corollary ID.21 and 
equation flD.ip together imply the following useful result. 

Corollary D.3. //Gp acts locally freely on then p is a regular value of the eguivariant 

momentum map fj. Conseguently, is an embedded submanifold of P. 


We hnish by a clarifying remark concerning the momentum map. 

Remark D.4. Let B be the Poisson tensor on P, for which {f,h} = B{df,dh) = Cv^f- 
We can write V/ = B^{dh) with the corresponding bundle map B'^\ T*P -E TP. Any 
X E Q = TeG = {Te'G*)* extends to a unique right-invariant 1-form {tx on G* {e E G and 
e! G G* are the unit elements). With this at hand, equation fID.ip can be reformulated as 


Xp = 


(D.5) 


which is a slight variation of the dehning equation of the momentum map found in |18] . 
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